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Abstract. — In this paper, we study a class of fractional optimal control problems. A necessary 
condition for the existence of an optimal control is provided in [2l [3l I23L 1241 134] and references 

f"~~. ■ therein. It is commonly given as the existence of a solution of a fractional Pontryagin's system 

and the proof is based on the introduction of a Lagrange multiplier. Assuming an additional 
condition on these problems, we suggest a new presentation of this result with a proof using only 
classical mathematical tools adapted to the fractional case: calculus of variations, Gronwall's 

C_J ' Lemma, Cauchy-Lipschitz Theorem, and stability under perturbations of differential equations. 

_J , In this paper, we furthermore provide a way in order to transit from a classical optimal 

"^j ' control problem to its fractional version via the Stanislavsky's formalism. We also solve a strict 

Cu I fractional example allowing to test numerical schemes. Finally, we state a fractional Noether's 

theorem giving the existence of an explicit constant of motion for fractional Pontryagin's systems 
admitting a symmetry. 

^ I Keywords: Optimal control; fractional calculus; Noether's theorem. 

(N ■ AMS Classification: 26A33; 49J15. 

Introduction 

m 

^^ ■ The control theory is the analysis of controlled dynamical systems. These systems are 

^N I varied: they can be differential, stochastic or discrete. The optimal control theory concerns 

the determination of a control optimizing a certain cost. Consequently, this theory is strongly 
connected to the 18* century classical mechanic (variational principles, Euler-Lagrange equa- 
S>^ I tions, etc., see [5], 126^ [38]). Since the second world war, this theory has a considerable devel- 

^ • opment and one can find applications in many domains: celestial mechanic 0, biology |12j . 

hydroelectricity |13j . economy [16^ I18|, I21j . etc. The subject is widely treated and one can 
find a lot of references on the subject, see for example |1H, I19|, \30\ I45j . 

The fractional calculus, i.e. the mathematical field dealing with the generalization of the 
derivative to any real order, plays an increasing role in many varied domains as economy |14j 
or probability |37|, I44j . Fractional derivatives also appear in many fields of Physics (see |31| ) : 
wave mechanic [4], viscoelasticity [6], thermodynamics [32j . fluid mechanic in heterogeneous 
media |27|, I47|, [48], etc. A natural question then arises: can we develop optimal control 
theories for fractional differential systems ? 

Recently, a subtopic of the fractional calculus gains importance: it concerns the variational 
principles on functionals involving fractional derivatives. This leads to the statement of 
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fractional Euler-Lagrange equations, see [H [Tl I42j . A direct consequence is the emer- 
gence of works concerning a particular class of fractional optimal control problems, see 
[21 [31 123|, I24|, I34j and references therein. These studies usually use the Lagrange multiplier 
technique allowing to write these problems as problems of optimization without constraint 
of augmented junctionals. With a calculus of variations, authors then obtain a necessary 
condition for the existence of an optimal control. This condition is commonly given as 
the existence of a solution of a system of fractional differential equations called fractional 
Pontryagin's system. 

In this paper, we first give a new presentation of this result. Precisely, making an additional 



assumption (see Condition (fx lip ) ) , we rewrite directly these fractional optimal control prob- 
lems as simpler problems of optimization without constraint of functionals depending only on 
the control. With a calculus of variations, we finally conclude to the same necessary condition 
but with a new presentation. Although the method here used is considerably inspired by 
the Lagrange multiplier technique, it allows us to give a complete proof using only classical 
mathematical tools adapted to the fractional case: calculus of variations, Gronwall's Lemma, 
Cauchy-Lipschitz Theorem and stability under perturbations of differential equations. 

Nevertheless, the explicit computation of controls satisfying the above necessary condition 
needs the resolution of the fractional Pontryagin's system which is a main drawback. Indeed, 
solving a fractional differential equation is in general very difficult. Consequently, in this 
paper, we suggest a deviously way in order to get informations on the solutions of a frac- 
tional Pontryagin's system. Precisely, we study the existence of classical conservation laws, 
i.e. functions which are constant on each solution. Indeed, constants of motion, generally 
associated to physical quantities, give strong informations on the solutions in the phase space 
for example. Moreover, they also can be used in order to reduce or integrate the system by 
quadrature. Previous results in this direction have been obtained by Torres and Frederico in 
|23|, I24j . However, in each of these papers, the conservation law is not explicit but implicitly 
defined by a functional relation. In this paper, inspired by a recent result obtained in |10j . 
we prove a fractional Noether's theorem providing an explicit conservation law for fractional 
Pontryagin's systems exhibiting a symmetry. 

In this paper, we also suggest: 

— a link between classical optimal control problems and their fractional versions via the 
Stanislavsky's formalism (see |33|, I44j ): 

— a solved fractional example allowing to test numerical schemes in the strict fractional 
case. 

The paper is organized as follows. Firstly, we remind some classical definitions and results 
concerning fractional calculus (Section [TJ. Then, we present the class of fractional optimal 
control problems studied and the usual strategy using the Lagrange multiplier (Section 12. ip . 
Then, assuming an additional condition on the problem, we develop a new proof leading to the 
same result with a different presentation, see Section 12.21 We then suggest a way in order to 
transit from a classical optimal control problem to its fractional version via the Stanislavsky's 
formalism (Section [3]). Subsequently, we suggest some examples with, in particular, a solved 
fractional example allowing to test numerical schemes in the strict fractional case (Section 
[3]). Finally, in Section [5l we state a fractional Noether's theorem for fractional Pontryagin's 
systems admitting a symmetry. Technical proofs of Lemmas are provided in Appendix 1X1 



1. Usual definitions and results about fractional calculus 

Let us introduce the following notations available in the whole paper. Let a < 6 be two reals, 
let d, m £ 'N* denote two dimensions and let || • || be the euclidean norm of M and M™. 

1.1. Fractional operators of Riemann-Liouville and Caputo. — The fractional 
calculus concerns the extension of the usual notion of derivative from non-negative integer 
orders to any real order. Since 1695, numerous notions of fractional derivatives emerge over 
the year, see |35|, I40|, I43j . In this paper, we only use the notions of fractional integrals and 
derivatives in the sense of Riemann-Liouville (1847) and Caputo (1967) whose definitions are 
recalled in this section. We refer to [35|, I40|, I43j for more details. 

Let g € 'i^'^([a, 6],M ) and a > 0. The left (resp. right) fractional integral in the sense of 
Riemann-Liouville with inferior limit a (resp. superior limit b) of order a of ^f is defined by: 

Vt G]a, b], I^g{t) :=-^ f (t- yT-^giy) dy (1) 

r(a) Ja 

respectively: 

Vi G [a, 6[, I1g{t) := -1- / (y - tr-^g{y) dy, (2) 

r(a) Jt 

where F denotes the Euler's Gamma function. Let us remind that I'^g (resp. I+g) is 
continuous and can be continuously extended by in t = a (resp. t = b). Let us note that 
I^g (resp. —I}^.g) coincides with the anti-derivative of g vanishing in t = a (resp. t = b). For 
a = 0, let Pg = 1° 5 = 9- 

Now, let us consider < a < 1. The left (resp. right) fractional derivative in the sense of 
Riemann-Liouville with inferior limit a (resp. superior limit b) of order a of 5 is then given 
by: 

Vt g]g,6], D'^git) := j^[l'S''g){t) (resp. Vt G [a, 6[, Dlg{t) := -^^[ll-g){t)), (3) 
provided that the right side terms are defined. 

In the Riemann-Liouville sense, the strict fractional derivative of a constant is not zero. 
Caputo then suggests the following definition. For < a < 1, the left (resp. right) fractional 
derivative in the sense of Caputo with inferior limit a (resp. superior limit b) of order a oi g 
is given by: 

VtG]a,6], ,D'lg{t) := DZ{g - g{a)){t) (resp. Vi G [a, 6[, ,D^g{t) := D'^{g - g{b)){t)' 

provided that the right side terms are defined. Let us note that if g{a) = (resp. g{b) = 0), 
then cD^g = D°Lg (resp. cD'^g = D^g). 

In the classical case a = 1, the fractional derivatives of Riemann-Liouville and Caputo both 
coincide with the classical derivative. Precisely, modulo a (—1) term in the right case, we 
have Dl = ^Dl = -D\ = -^D^ = d/dt. 
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Finally, let us remind the following important result. If g G ^^([a, 5],R ) and < a < 1, 
then we have: 

Vt G]a,6], cD^g{t) = I^T'^m and Vt G [a,b[, cD^git) = -/|-°<7(i), (5) 

where g denotes the derivative of g. Consequently, in this case, we have cD^g (resp. cD'^g) 
is continuous and can be continuously extended in i = a (resp. t = b). 

1.2. Some properties and results about the fractional operators. — Let us precise 
that the properties and results developed in this section are well-known and already proved 
in many references, see |35|, I43j for example. However, these results are not always exactly 
presented as we need. In this case, we give a detailed proof for the reader's convenience. 

First, we remind two basic results concerning the fractional integrals. They are proved in 
|35|, I43j both using the Fubini's theorem. The first one yields the semi-group property of the 
fractional integral operators: 

Property 1. — Let g G '^'^([a, 6],R'^) and a\, a2 > 0. Then, we have: 

Vt G [a, b], rj o rjg{t) = rj+^-g{t) and I^^ o /-^g(t) = /^^+"^5(i). (6) 

The following second property is occasionally called fractional integration by parts. It is very 
useful for calculus of variations involving fractional derivatives: 

Property 2. — For any gi, g2 G '^^{[a,b],W^) and any a > 0, we have: 

f iZgi ■g2dt= f gi- 7^52 dt. (7) 

J a J a 

Let US introduce the following convention: a function defined on ]a, b] (resp. [a, b[) is said to 
be an element of ^^{[a, 6], M"^) if and only if it is continuous on ]a, b] (resp. [a, b[) and if it can 
be continuously extended in t = a (resp. t = b). From Section Fl.il we can give the following 
examples: 

- for g G 'r°([a,6],M'^) and a > 0, we have I^g, I^g G ^°([a,6],M'^); 

- for 5 G 'r^([a,6],R'^) and < a < 1, we have ^D'^g, cD%g G '^^{[a,b\,W^). 

Now, we prove some results of composition between the left fractional operators. One can 
easily derive the analogous versions for the right ones. 

Property 3. — Let g G '^'^([a, 6],M'^) and < q < 1. Then, cDI o L'^g is an element of 
'i^^{[a,b],W^) and for any t G [a,b], we have: 

,D'^oI^g{t)=g{t). (8) 

Let us assume additionally that cD'^g G '^'^([a, 6],M'^). Then, /° o ^D"g is an element of 
^^{[a,b],W^) and for any t G [a,b], we have: 

L^o,Dlg{t)=g{t)-g{a). (9) 

Proof. — Let us prove the first result. Since L"g{a) = 0, we have for any t G]a, b]: 

,D1 o I-g{t) = D^_ o I'^g{t) = ^ o /I- o L'^g{t) = j^[l-9){t) = g{t). (10) 

Hence, c-D" o I'^g is continuous on ]a, b] and can be continuously extended by g{a) in i = a. 
Now, let us prove the second result. It is obvious for a = 1. Now, let us consider < a < 1. 



Since ^D'^g e '^"([a, 6],]R'^), we have li""(ff - g{a)) G '^^([a, 6],R'^). Combining the first 
result and Equality ([5]), we have for any t €]a,b]: 

git) -g{a) = ,D\-'- o J^" [g _ ^(a)) (i) = r o A o /l"" (5 - 5(a)) (t) 

= roD'l{g- g{a)) (t) = II o ,D^^g{t). (11) 

Hence /° o cD^g is continuous on ]o, 6] and can be continuously extended by in t = a. The 
proof is completed. D 

Finally, we prove the following fractional Cauchy-Lipschitz type theorem. Let us precise that 
a version of this theorem is proved in |35l Part 3.5.1, Corollary 3.26, p. 205]. It gives the 
existence and the uniqueness of a global solution of a fractional differential equation. 

Theorem 4 (Fractional Cauchy-Lipschitz theorem). — Let < a < 1, F G 'tf^{M.'^ x 
[a,6],M ) and A G M . Let us assume that F satisfies the following Lipschitz type condition: 

3K £M, V(xi,X2,t) G (R'^f X [a,b], \\F{xi,t) - F{x2,t)\\ <K\\xi-X2\\. (12) 

Then, the following fractional Cauchy problem: 

,D'^g = F{g,t) 

g{a) = A, ^^^^ 

has an unique solution in 'ig^{[a, 6],M ). The solution g is an element of'rf^'^'{[a, b],M. ) (where 
[a] denotes the floor of a) and g satisfies: 

VtG[a,6], g{t) = A + I^{F{g,t))it). (14) 

Proof. — Firstly, let g G ^"^([a, 6],M ). From Property [3l we conclude that g is solution of 
the fractional Cauchy problem ()13p if and only if g satisfies Equality ()14p . Now, let us define: 

f. '^°([a,6],]R"') -^ ^°([a,6],R'^) (15) 

g ^ ^{g): [a,b] -^ R'^ 

t ^ A + I^{F{g,t)){t). 

By induction, we prove that for any n G N* and any gi, g2 ^ '^^^cl-, ^]>]R'^)i we have: 

Vi G [a, 6], \W^{g{)(t) - (^"(<72)(t)|| < T(TT^"^' " ^'"°°- ^^^^ 

Consequently, we have for any n G N* and any 51, 52 £ '^'^([a, 6],M''): 

w{9i) - 9^"(52)iioo < ^;" \ ibi - 52II00. (17) 

r(l + no) 

Since K'^(\) — a)"'"/r(l + no) — > 0, there exists n G N* such that (/?"■ is a contraction. 
Consequently, 93 admits an unique fix point in "^^([a, 6],M'^). The proof is completed. D 



2. A class of fractional optimal control problems 

From now and for all the rest of the paper, we consider < q < 1 and A G 
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2.1. Presentation of the problem and usual strategy. — In this section, we first give 
a brief presentation of the class of fractional optimal control problems interesting us in this 
paper. Finally, we briefly remind the usual strategy (using the Lagrange multiplier technique) 
leading to a necessary condition for the existence of an optimal control, see [2], [3l I23|, I24|, I34j 
and references therein. 

An optimal control problem concerns the optimization of a quantity depending on parameters 
given by a controlled system. Precisely, the aim is to find a control optimizing a certain cost. 
Such a control is called optimal control. Many of these studies lead to the use of the Lagrange 
multiplier technique. For example, we refer to (SI I20|, I28|, 129] for controlled systems coded 
by ordinary differential equations. 

In this paper, we work in the following framework. We are interested in systems controlled 
by the following fractional Cauchy problem: 

,D":q = f{q,u,t) , . 

q{a)=A, ^ > 

where u denotes the control. Our aim is to find a control u optimizing a quantity of the form: 

eh 

L{q,u,t)dt, (19) 



where q is solution of (fT8|) . 

The common strategy is first to write this problem as a problem of optimization under con- 
straint: 

nb 

argmin / L{q,u,t) dt. (20) 

{q,u) satisfying I I18II J a 

Then, the Lagrange multiplier technique consists in the study of the critical points of the 
following augmented functional: 

l-b 
{q, u,p)^^ / L{q, u,t)-p- {cDlq - f{q, u, t)) dt, (21) 

J a 

where p is commonly called Lagrange multiplier. Finally, with a calculus of variations, such a 
strategy leads to the following result: a necessary condition for {q, u) to be a solution of ()20p 
is that there exists a function p such that the following fractional Pontryagin's system holds: 

( dH 

cD'^q = —{q,u,p,t) 

D%P=^{q,u,p,t) ^^^^ 

— {q,u,p,t) = 

{q{a),p{b)) = {A,0), 

where H{x, v, w, t) = L{x, v,t) +w ■ f{x, v, t). We refer to [2|, [3|, I23|, I24^ I34j and references 
therein for more details. 



In this paper, we are going to assume that / satisfies a Lipschitz type condition, see Condition 



( fx hp ) . This assumption will allow us to write directly the initial problem as a simpler prob- 
lem of optimization without constraint of the initial functional which is then only dependent 
of the control u. Finally, a simple calculus of variations leads us to the same result but with a 
new presentation. Let us remind that the method here developed is widely inspired from the 
Lagrange multiplier technique. Nevertheless, it allows us to give a complete proof only using 
classical mathematical tools adapted to the fractional case: calculus of variations, GronwaU's 
Lemma, Cauchy-Lipschitz Theorem and stability under perturbations of differential equations. 

2.2. New presentation of the result. — In this section, we first give rigorously the 
definitions concerning the class of fractional optimal control problems described in Section 

ED 

— The elements denoted u G 'rf^{[a,b],W"^) are called controls; 

— Let / be a '^^ function of the form: 

/ : M'^ X M"' X [a, b] — > M"^ (23) 

{x,v,t) I — > f{x,v,t). 

It is commonly called the constraint function. We assume that / satisfies the following 
Lipschitz type condition. There exists M > such that: 

y{xi,X2,v,t) G (R'^)2 X M"* X [a,b], \\f{xi,v,t)-fix2,v,t)\\ < M||xi -xall; (/. lip) 

— For any control u, let g"'" G '(^["^([a, 6],M'^) denote the unique global solution of the 
following fractional Cauchy problem: 

' cD'^q = f{q,u,t) 



qia) = A. (^^9 ) 

qU,a jg commonly called the state variable associated to u. Its existence and its unique- 



ness are provided by Theorem U] and Condition (fx lip); 
— Finally, with this condition on /, the fractional optimal control problem described in 
Section 12.11 can be rewritten as the simpler problem of optimization of the following cost 
functional: 

C" : ^0([a,6],M"^) — > R (24) 

u ^^ I L{q''''^,u,t)dt, 

J a 

where L is a Lagrangian, i.e. a ^^ application of the form: 

L: M°' X M™ X [a, 6] — > M (25) 

{x,v,t) I — > L{x,v,t). 

Hence, the existence and the uniqueness of q"'° for any control u allow us to rewrite directly 
the initial problem as a simpler problem of optimization without constraint of the cost 
functional £°: we do not need to introduce an augmented functional with a Lagrange 
multiplier. Moreover, let us note that C°' is only dependent of the control. 

A control optimizing £" is called optimal control. A necessary condition for a control u to be 
optimal is to be a critical point of £°, i.e. to satisfy: 

VnG^O([a,&],M"*), Z)£-(n)(n) := lim ^^^^^^^^^^^^— ^^ = 0. (26) 

£-s>o e 
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In the following, we then focus on the characterization of the critical points of >C". With an 
usual calculus of variations, we obtain the following Lemma [5] giving explicitly the value of 
the Gateaux derivative of £": 

Lemma 5. — Let u, u £ 'rf^{[a,b],W^). Then, the following equality holds: 

— ((?"'",«, t).q+ — (g"'", u,t).u dt, (27) 

where q G '^^"'([a, 6],IR ) is the unique global solution of the following linearised Cauchy 
problem: 



c 



df df 

D''q= -^(g"'",n,t) xg+-^(g"'°,n,t) xu 
"^ dx^^ ' ' ^ ^ dv^^ ' ' ' (LCPf) 

/(«) = 0. 
Proof. — See Appendix IA.4I D 

This last result does not lead to a characterization of the critical points of £" yet. Then, let 
us introduce the following elements stemming from the Lagrange multiplier technique: 

— Let H he the following application 

H : M'' X M™ X M'^ X [a, 6] — > R (28) 

{x,v,w,t) I — > L{x,v,t) + w ■ f{x,v,t). 

H is commonly called the Hamiltonian associated to the Lagrangian L and the constraint 
function /; 

— For any control u, let p"'" G '(^["l([a, ^JjIR"^) denote the unique global solution of the 
following fractional Cauchy problem: 

ji(6) = 0. 

pU,a jg commonly called the adjoint variable associated to u. Its existence and its 
uniqueness are provided by the analogous version of Theorem U] for right fractional 
derivative. Since p"'"(6) = 0, we can write cD^p^'" = D'^p'^'"'. 

Let us note that, for any control u, the couple (g^'^^p"'") is solution of the following fractional 
Hamiltonian system: 

cD'^q = -^-{q,u,p,t) 

D'^p= —{q,u,p,t). 
Finally, the introduction of these last elements allows us to prove the following theorem: 

Theorem 6. — Let u G 'tf^{[a,b],W^). Then, u is a critical point of CJ^ if and only if 
(g"'", u,p"'") is solution of the following fractional stationary equation: 

f)J-f 

— {q,u,p,t)=0. (SE") 



Proof. — Let u, u ^ 'i^^{[a,b],W^). From LemmaEl we have: 

— {q-'-,u,t)-q+—{q-'^,u,t)-udt. 



(29) 



Then: 



°'^°*"'<"' = r (i<'"'°'"'" + (!<'"'"• "•'))' " '"") 



- ( If (9"'°.".') X f) ■!'°" + f;(«"'°.".') ■ » <«*■ (30) 



From Theorem m since q is solution of (LCP^^) and p"'"^ is solution of (CP" ), we have: 



I-{^{q-'-,u,t)xq+^{q^'^,u,t)xn 



and 



P 



u,a T(x 



^+ (§('^"'"'"'')+ (If ^'^"'"'"''0^ ''''"'") 



(31) 



(32) 



Then, using the fractional integration by parts given in Property O we obtain: 
rb 



DC''{u){u) 



P 



^{q-^-,u,t)xq + ^{q^''^,u,t)xu 



Then: 



Z)/:"(n)(n) 



^{q-^'-,u,t)xq\ .p-^- + ^{q-^-,u,t)-udt. (33) 



^iq-'-,u,t)] xp-'" + |^(g"'",u,t)] -udt 



dH , 

dv 



{q''''',u,p''''',t)-udt. 



The proof is completed by the Dubois-Raymond's lemma. 



D 



Let us note that Lemma [5] is proved in Appendix lA. 41 using Lemmas ll6lll7l andll8l These last 
three Lemmas are respectively a fractional Gronwall's Lemma, a result of stability of order 



1 and a result of stability of order 2 for the fractional Cauchy problem (CP" ). Hence, the 
proof of Theorem[6]is only based on classical mathematical tools adapted to the fractional case. 

From Theorem [6l we retrieve the following result provided in [2], \3\. \23[ I24|, [34] and references 
therein: 



Corollary 7. — £" has a critical point in^^{[a,b],W^) if and only if there exists {q,u,p) £ 
^["] ([a, b] , R'^) X '^°([a, 6] , M™) x '^["1 ([a, b] , W^) solution of the following fractional Pontryagin's 
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system: 

BJ-f 

cD'^q= —{q,u,p,t) 

BTJ 

D%p= —{q,u,p,t) 

dH , , ^ 

— {q,u,p,t) =0 

[ {q{a),p{b))={A,0). 
In the affirmative case, u is a critical point of C" and we have {q,p) = (g"'",p"'") 



(PS°) 



In practice, we are going to use more Corollary [7] than Theorem [6l see Examples in Section 
m Let us note that the fractional Pontryagin's system ()PS"P is made up of the fractional 
Hamiltonian system (|HS"|) . the fractional stationary equation (|SE"P and initial and final 
conditions. 



2.3. Remark 1: the classical case. — In the case a = 1, the fractional derivatives 
coincide with the classical one. Consequently, the fractional optimal control problem studied 
coincides with the classical one. Then, in this case, Corollary [7] is nothing else but the classical 
theorem obtained in P [20l [28l [29] : 

Theorem 8. — C^ has a critical point in'^^{[a,h],W^) if and only if there exists {q,u,p) £ 
^i([a,6],M'^) X ^0([a,6],M™) x '^^{[a,b],W^) solution of the following Pontryagin's system: 

9 = ^r-W'^'^'O 



dw 
dH 

dx 



{q,u,p,t) 



dH 

dv 



(PS'] 



{q,u,p,t) = 



I {q{a),p{b))={A,0). 
In the affirmative case, u is a critical point of C^ and we have {q,p) = {q^'^,p^'^). 



2.4. Remark 2: link with the fractional Euler-Lagrange equation. — Let us take 



the constraint function f{x,v,t) = v satisfying (fx lip). In this case, the fractional optimal 
control problem studied is: 



optimize / L{q, u, t) dt under the constraint <j '^ / T _ /. 



(34) 



Using Corollary [71 if this problem has a solution, then there exists a solution (q,u,p) £ 
^[°l([a,6],M'^) X <r°([a,6],M'^) x '^^°'\[a,b],M.'^) of the fractional Pontryagin's system ([P^^ 
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here given by: 

( cD'^q = u 

dL 

Dlp=—{q,u,t) 

^^ (35) 

— {q,u,t) +p = 

( {q{a),p{b))={A,0). 

In the case of the existence of an optimal control u for problem p4p , we obtain that the state 
variable associated g"'" is solution of the following fractional Euler- Lagrange equation: 



— {q,,D''_q,t)+D'l^—{q,,D''_q,t)j=0. (EL^ 



According to the works of Agrawal in [T], we then obtain that q^'" is a critical point of the 
following fractional Lagrangian functional: 



L{q,^Dtq,t)dt. (36) 

We refer to [1] for more details concerning fractional Euler-Lagrange equations. 

3. A transition from the classical to the fractional problem 

As seen in Introduction, there is a large development in fractional calculus and consequently 
concerning optimal control problems with controlled systems coded by fractional differential 
equations. In this section, we present a link between a classical optimal control problem and 
its fractional version via the Stanislavsky's formalism, see |33|, I44j . Only for this section, we 
assume that [a, h\ = [0, r] with r > and < a < 1. 

3.1. The Stanislavsky's formalism. — The governing equation for a fluid in a ho- 
mogeneous porous medium is usually the classical Richards equation and it is derived 
from physical principles, see |41j . Nevertheless, the fractional Richards equation plays 
an important role in the study of the behaviour of a fluid in a heterogeneous medium, 
|27|, I47^ I48| . A natural question then arises: how physically understand the emergence of 
the fractional order when the medium changes? Let us sketch an example of answer. A 
heterogeneous porous medium presents a lot of stochastic heterogeneities implying complex 
geometric structure and then inducing stochastic retention zones, |47j . Hence, we can 
assume that a change of medium induces a modification of the time variable. Therefore, 
following the work of Stanislavsky consisting in the introduction of a stochastic internal time 
(a "slow" time), we remind that the emergence of a fractional derivative can be the effect of 
a change of time. Let us give more details. However, we refer to |33|, I44j for a complete study. 

Let us consider Tt a stochastic process of probability density function p{y, t). Tt represents the 
stochastic internal time (the "slow" time). Stanislavsky assumes that the Laplace transform 
Lap(„) of /?(y, t) with respect to its first variable satisfies: 

Vt > 0, Vs G R, Lap(j^) [p(y , t)] (s) = E«,i (-si") , (37) 
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where E^^i is the Mittag-Leffler function defined in Appendix lA.li A possible construction 
of Tt is given in 



Stanislavsky studies the dynamical effects of this change of time. Precisely, for any g G 
^'^([0,r],M'^) and any t G [0,r], he studies the function Ta{g) defined by: 

TU9m=^g{Tt)), (38) 

where E designates the mean value. The main property is the following result proved in 



Lemma 9. — Let g ^ '^^([0,r],M'^). Then, we have: 

Vt G [0,r], ^g{Tt)) = ,Dl{F^{g)){t). (39) 

Hence, a classical derivative is transformed into a fractional one under the introduction of the 
stochastic internal time of Stanislavsky. 

3.2. Application on the class of fractional optimal control problems. — Now, let 
us talk about the consequences of the works of Stanislavsky on the class of fractional optimal 
control problems studied in this paper. We use the notations and definitions of Section [2.21 



Let us assume that L and / are autonomous, i.e. L{x,v,t) = L{x,v) and f{x,v,t) = f{x,v). 
Then, let us assume that / satisfies the following condition: 

y{q,u) G '^\[0,t],R'^) X ^°([0,r],]R'"), Vt G [0,t], 

E(^f{q{Tt),uiTt))) = f(E{qiTt)),E{u{Tt))). (40) 

Such a condition is satisfied by linear autonomous constraint functions, as the examples in 
Sectional According to Lemma [9l one can prove that for any control u, Fa{q^'^) satisfies: 

r c^^j-„(g-'i) = /(j-,((?"'i), j-,(n)) 

\ To.{q^^^){a) = A. ^^') 

Precisely, we have for any control u, J^a{q^'^) = q-^"('^)'". 

Now, let us consider the optimal control problem studied in Section 12.21 in the classical case: 

optimize the cost functional u G ^°([a, &],]R"') i — > L{q'''^,u,t) dt. (42) 

Jo 

Let us assume that the controlled system can be prone to changes inducing a modification of 

the time variable as it is the case with porous media and the Richards equation. We are then 

interested in the optimization of the following "slow" cost functional: 

nG'^°([o,6],M^)^ / L(j«(g"'i),J-„(M)) dt (43) 

Jo 

which is relied to the following fractional cost functional: 

wG'^°([a,6],]R"')^ rL(g"'",n)dt (44) 

Jo 



as its restriction to the image of J-q,. Hence, once an optimal control of the functional 
found, one can be interested in its projection on the image of J-q, in order to approach an 
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optimal control of the functional (03]) . 

Hence, the Stanislavsky's formalism is an example relying a classical optimal control problem 
to its fractional version via the introduction of a stochastic internal time. 



4. Examples 

In this section, we are going to study some examples of optimal control problems studied in 
Section [2.21 both in classical and fractional cases. 



4.1. An applied classical linear-quadratic example. — Classical linear-quadratic 
examples are often studied in the literature because they are used for tracking problems. The 
aim of these problems is to determine a control allowing to approach as much as possible 
reference trajectories, [451 Part 1.4, p. 49]. In this section, we study such an example, |191 
Part 4.4.3, example 3, p. 53]. More generally, a quadratic Lagrangian is often natural (for 
example in order to minimize distances) and even if the differential equations are frequently 
non linear, we are often leaded to study linearised versions. 

In this section, we take a = d = m = 1. We consider an evolution problem of two nondescript 
populations z and ^ interacting during a time interval [a,b]. We control the injection or the 
discharge of the population ^ in the system: we then control the value of the population ^ 
in real time. The interaction between z and ^ is governed by the following linear differential 
equation: 

z = z + ^. (45) 

For any control S, and any real zq, we denote by Zz^,^ the unique solution of (I45p associated 
to ^ and satisfying the initial condition z{a) = zq G M. Let us assume that, for an initial 
condition z{a) = X G M, we know experimentally a satisfactory control ^x such that the 
evolution of the population zx,^x ^^ healthy with respect to a nondescript constraint. 

The problem is finally the following: assume that the initial condition is modified (i.e. z{a) = 
Y 7^ X), what is the control ^ minimizing the difference between ^ and ^x plus the difference 
between zy,s, cLnd zx,^x ^ More precisely, we are looking for a minimizer of the following 
functional: 

e G ^\[a, 6],M) ^ 1 ^ (e - ix? + {zY^i - zx,u? dt. (46) 

With a change of variable A = Y — X^ u = ^ — ^x-, Q = zy^^ — zx,^x ^^^ giving the following 
quadratic Lagrangian and the following linear constraint function: 

L: M2 X [a, 6] — > R and / : M^ x [a, 6] — > R (47) 

{x,v,t) I — > (x^ + -y^)/2 {x,v,t) I — > x + v, 

the problem is reduced to find an optimal control u for the cost functional C^ associated to 
L, / and A. See Section 12.21 for notations and definitions. 
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According to Theorem [HI we are interested in solving the Pontryagin's system (PSi) here 
given by: 

q = q + u 

p = -q-p 

p + u = 
[ (g(a),p(6)) = (A0). 
Hence, with p + u = 0, it is sufficient to solve the following linear Cauchy problem: 



(48) 



q = q + u 

u = q- u (49) 

(g(a),n(6))=(A0). 

Using matrix exponentials, one can prove that it exists an unique solution (q, u) to Cauchy 
problem (|49p given by: 



VtG [a, 6], 



l-R 



q{t) = A cosh (\/2(t - a)) H -=- sinh [y/2{t - a)) 



u{t) = A 



\l + R 



V2 



^/2 



sinh (\/2(t - a)) - Rcosh {V2{t - a)) 



where 



it- 



sinh {y/2(b - a)) 



(50) 



(51) 



V2cosh {V2{b - a)) - sinh {^/2{b - a)) ' 

Finally, we conclude that u is the unique critical point of C^. Hence, if there exists S^ mini- 
mizing (|46p . then: 

yte[a,b], at)=^x(.t) + u{t). (52) 

4.2. The fractional linear quadratic example. — In this section, we consider the 
previous example of Section 14.11 in the strict fractional case < a < 1. Let us consider 
d = m = 1 and the couple (L, /) given in Equation (j47p . We want to express a critical point 
u of the cost functional C" associated to L and /. See Section [2^2] for notations and definitions. 

Therefore, according to Corollary [TJ we look for (q, u,p) solution of the fractional Pontryagin's 
system ()PS"P here given by: 

cD°_q = q + u 

D%p = q+p 

p + u = 
[ {q{a),p{b))={A,0). 

Hence, from n + p = 0, it is sufficient to solve the following linear fractional Cauchy problem: 

cDlq = q + u 

D%u = u-q (54) 

{q{a),u{b))=iA,0). 

Although the fractional Cauchy problem (|54p is linear, the presence of the left and the 
right fractional derivatives is a main drawback to the explicit computation of a solution. 
More generally, this characteristic implies very big difficulties in order to solve the most of 
fractional Hamiltonian systems. Finally, the resolution of (|PS"P is the most of time still an 
opened problem. 



(53) 
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Nevertheless, in this case, we can hope getting informations on the solutions of ()PS"P in 
dimension d = m = 2 with the help of the fractional Noether's theorem proved in Section 
[5l Indeed, in this case, the fractional Pontryagin's system (|PS"|) admits a symmetry and a 
constant of motion can be obtained. We refer to Section [5] for more details. 

4.3. A solved fractional example. — In this section, we give a solved fractional exam- 
ple in the sense that we give explicitly the unique critical point of the cost functional £". 
Nevertheless, the fractional Pontryagin's system ()PS°P is still not completely resolved. 

Let us take d = m = 1, [a,b] = [0, 1] and < a < 1. Let us consider the following Lagrangian 
L and the following linear constraint function /: 

L: M2 X [0, 1] — > R and / : M^ ^ [o, 1] — > R (55) 

ix,v,t) I — > {v'^/2) + -f{l-t)^x {x,v,t) I — > Xx + fiv, 

where /3, 7, //, A are parameters in M*. We are looking for a critical point u of C" associated 
to L and /. See Section [2.21 for notations and definitions. 

According to Corollary [71 we have to solve the fractional Pontryagin's system (|PS"P here 
given by: 

cD'^q = Xq + fiu 

u + fj.p = 

((?(o),p(i)) = (Ao). 

From Theorem U the Cauchy problem cD'^q = Xq + fiu with q{0) = A admits an unique 
solution for any u G '^'^([0, 1],M). Then, from u + fip = 0, it is sufficient to solve the following 
fractional Cauchy problem: 

DIP = Xp + 7(1 - tf ^q^i^^i^^t to I ^-P' = ^P^ + ^'^ (57) 

p{l) = equivalent to | ^^^^^ ^ ^^ [^i) 

with the change of unknown po{t) = p{l — t) for any t € [0, 1]. The unique solution of this 
last fractional Cauchy problem is given in |35[ Chap. 3, p. 137] by: 

Vt G [0, 1], po{t) = f {t- y)"-iE„,,(A(t - yr)jy^ dy, (58) 

Jo 

where J^a,a is the Mittag-Leffler function defined in Appendix lA.il In order to get a better 

formulation, we make a change of variable which gives us: 

Vt G [0, 1], po{t) = 7 f {t- yfy"-^E^,^{Xy'')dy (59) 



Jo 

= 7r(/3 + 1)/^' (y"-iE„,,(Ay")) (t) (60) 

= 7r(/3 + l)t"+^E„,,+/3+i(At'^). (61) 

We refer to |35j for more details concerning the calculations. Finally, we obtain the unique 
critical point u of £" given by: 

Vt G [0, 1] , u{t) = -fijrifS + 1) (1 - tr+^E^,a+^+i ( A(l - t)") . (62) 

Although we obtain the unique critical point u of £", let us note that this example does 
not provide a completely solved fractional Pontryagin's system: the state variable g"'" is 
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still unknown. However, this example allows to test the quality of numerical schemes giving 
approximations of critical points of cost functionals C^ . This will be done in a further paper 
for a variational integrator of (|PS"p . 

5. A fractional Noether's theorem 

As we have seen with concrete examples in Sections 14.21 and 14.31 fractional Pontryagin's 
systems (jPS"P are the most of time not resolvable explicitly. This is a strong obstruction in 
order to express explicitly a critical point of a cost functional £°. In this section, we suggest 
a deviously way in order to get informations on the solutions of ()PS°P and consequently on 
the critical points of the cost functional £° associated. Precisely, we are interested in the 
existence of conservation laws for fractional Pontryagin's systems admitting a symmetry (see 
Definition dlj. 

Indeed, in 1918, Noether proved the existence of an explicit conservation law for any classical 
Euler-Lagrange equation admitting a symmetry. We refer to [5], I36^ 139] for more details. 
Adapting her strategy to the fractional case, Cresson, Torres and Frederico proved in 2007 
a preliminary result giving a conservation law for any fractional Euler-Lagrange equation 
admitting a symmetry, see |15^ I22^ I25| . Nevertheless, the conservation law obtained was 
not explicit but only given implicitly via a functional relation. From this first result, we have 
formulated in |10] a fractional Noether's theorem providing an explicit formulation of this 
conservation law via a transfer formula. 



At the same time, in |23|, [24], Torres and Frederico applied a similar strategy for fractional 
Pontryagin's system admitting a symmetry. Nevertheless, this conservation law is also given 
iinplicitly via a functional relation. In this section, we are going to formulate a fractional 
Noether's theorem providing an explicit formulation of this conservation law via an other 
transfer formula. 

We first review the definition of a one parameter group of diffeomorphisms: 

Definition 10. — Let n G N*. For any real s, let 4>{s,-) : M" — > M" be a diffeo- 
morphism. Then, $ = {(/>(s, •)}sgiR is a one parameter group of diffeomorphisms o/M" if it 
satisfies: 

1. ^{0,-) = IdK^; 

2. Vs, s' G M, ^{s, •) o (f){s', •) = ^{s + s', ■); 

3. (p ^-5 of class '^'^ . 

Usual examples of one paraineter groups of diffeomorphisms are given by translations and 
rotations. The action of three one parameter groups of diffeomorphisms on an Hamiltonian 
allows to define the notion of a symmetry for a fractional Pontryagin's system: 

Definition 11. — For i = 1,2,3, let <!>,, = {0i(s, •)}sgK be a one parameter group of diffeo- 
morphisms o/M"^, M™ andM.'^ respectively. Let L be a Lagrangian, f be a constraint function 
and H he the associated Hamiltonian. H is said to he cD'^-invariant under the action of 
{^i)i=i^2,3 if it satisfies for any {q,u,p) solution of (jPS°P and any s £ M.' 

H[4'i{s,q),4>2{s,u),(i)3{s,p),tj - 03(5, P) •c-D"((/>i(s,g)) = H{q,u,p,t) -p-cD'lq. (63) 
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From this notion, Torres and Prederico proved in 



the following result: 



Lemma 12. — Let L he a Lagrangian, f be a constraint function and H be the associated 
Hamiltonian. Let us assume that H is ^D'^-invariant under the action of three one parameter 
groups of diffeomorphisms {^i)i=i2 3- Then, the following equality holds for any solution 

iq,n,p) of mS^ 

d(t>i 



■Dt 



(t(»-))--t(»-)->-' 



(64) 



Then, from this result, Torres and Frederico defined a notion of fractional- conserved quantity. 
Nevertheless, this result did not provide exactly a constant of motion. In this section, using a 
transfer formula, we are going to write the left term of Equation (j64|) as an explicit classical 
derivative and then we obtain a real constant of motion. Let us provide this transfer formula: 



Lemma 13 (Transfer formula). — Let gi, g2 G 'tf°°{[a,b],W^) satisfying the following 
condition (C); 

the sequences of functions (/^~"' 



M\ 



Ap) 



[91- 9i{a)) -sDpgN* "-^^ 



{gf' • /^ '^5'2)pGN* converge uniformly to on [a,b\. 
Then, the following equality holds: 

d 



.D'lgi ■92-91- Dlg2 



dt 



Y^i-iyr 



r+l—a 



91 



W , Ar) jr+l-a 



9i{a))-9r+9r-r 



92 



.r=0 



Proof — In [lOj . under a similar condition than Condition (C), we have proved: 



^-ffi • 92 



d 
di 



Y^{-iyr_+^-g,.g^^ 



(r) 



,r=0 



and - gi ■ £'+52 



d_ 
It 



Hi 

,r=0 



(r) jr+l-a 
■ -'+ 92 



(65) 



(66) 



Consequently, writing cT>°Lgi = T)°L[gi — gi{a)), the proof is completed. D 

A discussion is provided in |10j concerning the condition (C): one can prove that this condition 
is satisfied for any couple of analytic functions for example. Then, combining Lemmas 1121 and 
[THl we prove: 

Theorem I4 (A fractional Noether's theorem). — Let L be a Lagrangian, f be a con- 
straint function and H be the associated Hamiltonian. Let us assume that H is ^D'^-invariant 
under the action of three one parameter groups of diffeomorphisms (*I'j)i=i,2,3- Let {q,u,p) be 
a solution of ()PS"P and let g denote d(j)i/ds{0,q). If g and p satisfy Condition (C), then the 
following equality holds: 



d 
di 



^(_l)r-r+l-a (^ _ ^(^)) . pir) ^ gir) . jr+l-y 



r=0 



0. 



(67) 



This theorem then provides an explicit algorithmic way to compute a constant of motion for 
any fractional Pontryagin's systems admitting a symmetry. Nevertheless, the conservation 
law is only given as a series of functions: in most cases, it is not easily computable. However, 
an arbitrary closed approximation of this quantity can be obtained with a truncation. 

Let us note that the fractional linear-quadratic example developed in Section 14.21 gives a 
concrete example of fractional Pontryagin's system admitting a symmetry: 
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Example 15. — Let us consider d = m = 2, the following quadratic Lagrangian and the 
following linear constraint function 

L: M2 X m2 ^ [0, 1] — > R and f : M^ x M^ x [0, 1] — > M^ 

{x,v,t) I — > (||x|p + ||ti|p)/2 {x,v,t) I — > x + v. 

(68) 
Then, let us consider the three one parameter groups of diffeomorphisms given by the following 
rotations: 

X IR2 ^ ^2 (gg) 



(. r, ro^ ^ ( cos(s^./) -sm{s9i) \ ( xi\ 

{S,Xl,X2) ^ ^^ sm{se,) COsisOi) J \ X2 J 



for i = 1, 2, 3 and where 6i, 02 G I^ cmd 9^ = —9i. With these parameters, one can prove that 
the Hamiltonian H associated to L and f is ^Dl-invariant under the action of (^j)j=i,2,3- 
Consequently, the fractional Pontryagin's system given in ()53p (in dimension 2) is not re- 
solvable but admits a symmetry and then admits an explicit conservation law given by the 



fractional Noether's Theorem 14' As said previously, this constant of motion is not explic- 
itly computable. However, truncating the series, one can provide an approximation of this 
quantity. 

A 
Appendix 

A.l. A fractional Gronwall's lemma. — Let us recall the definition extracted from 
|35|, I40|, I43| of the Mittag-Leffler function Eq,^ 0,3 of parameters ai, 02 > 0: 

00 f^ 

Vt G M, E,,,,,(t) := Y, r( u ^ V (^0) 

^^T{aik + a2) 

Now, let us give the following Lemma proved in |17|, I46j . For the reader's convenience, we 
recall the proof: 

Lemma 16 (A fractional Gronwall's lemma). — Let g G '^"([a, 6],]R) and a > 0. Let 

us assume that there exist K\, K2 > such that: 

Vt e [a, 6], < g{t) < Kil^g{t) + K2. (71) 

Then, g satisfies: 

yt£[a,b], 0<g{t)< K2Ea,i {Ki {t - a)") . (72) 

Proof. — Using Property [H Assumption ()7ip implies by induction that for any n G N* and 

any t £ [a, b]: 

n— 1 n— 1 / jy- (j. _ \a\^ 

< g{t) < K^rj-git) + K2Y, K^ll^iim = K^P''g{t) + K2 Y. viak + l) ' ^^^^ 

fe=0 A:=0 ^" ' 

Since g is continuous on [a, b] and non-negative, g is bounded by and by a constant K^ > 0. 
Hence, using the Stirling asymptotic formula, we have for any t G [a, b]: 

K"- f* (K^(t — a)"]^ 

< K^r^git) = — ^ / {t - yr~'g{y)dy < K, ^ '^ \^ -, 0. (74) 
1 [an) J^ 1 [an + Ij n^+00 

Finally, making n tend to 00 in Inequality ()73p . the proof is completed. D 
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Let us note that Ei^i is the exponential function. Consequently, in the classical case a = 1, 
the fractional Gronwall's lemma is nothing else but the classical Gronwall's lemma. 

A. 2. Result of stability of order 1. — In this section, we use the notations and definitions 
given in Section 12.21 We prove the following Lemma with the help of Lemma [T6l 



Lemma 17. — Let u, u ^ 'io^{[a,h],W"'). Then, there exists a constant Ci > such that: 
V|e| < 1, Vt G [a, 6], ||g"+""'"(t) - g"'"(t)|| < Ci\e\. (75) 

Proof. — From Theorem HI we have for any |e| < 1 and any t S [a, 6]: 

g"'°(t) = A + r {fiq'''",u, t)) (t) and g"+^"'°(t) = A + I^ (/(g"+^"'", u + eu, t)) (t). (76) 
Hence, for any |e| < 1 and any t E [a, b], we have: 



u+eu,a 



{t)-q-'^ 



l^a (/(g«+^^,«, u + eu, t) - fiq^'^, u, t)) (t) II 



< /- (11/(9"+^"'", « + eu, t) - f{q-'», u, t)\\ ) (t). 
We have for any |e| < 1 and any y G [a, h\. 



(77) 
(78) 



||/(g"+^^'"(y), u{y) + eu{y),y) - f{q^-''{y), u{y),y) \\ 

< \\f{q''^''''"{y),u{y)+eu{y),y)-f{q'^'^{y),u{y)+eu{y),y)\\ 

+ ll/('?"'"(y),n(y)+en(y),y)-/((?"'"(y),n(y),2/)||. (79) 



From Condition (f^ lip) and with a Taylor's expansion of order 1 with explicit remainder, the 
following inequality holds for any |e| < 1 and any y £ [a, b]: 



\\f{q-+^^^-iy),u{y)+euiy),y) - f{q-'''iy),u{y),y)\\ 
<M||g"+^"'"(y)-g"'"(y)|| + |e| 



|{(g"'"(y),f(y),y)xu(y) 



(80) 



where ^^(y) G [u{y),u{y) + £u{y)] C [—Mi, Mi]™ with Mi > independent of |e| < 1 since u 
and u are continuous on [a,b]. 

Thus, since df/dv, q'^'" and u are continuous, there exists M2 > such that for any |e| < 1, 
we have: 



Vy G [a,b], 



|{(g"'"(y),f(y),y)xu(y) 



<Mo 



|^(g"'",f,t)xn 



and then, Vt G [a,b], /' 
where M3 := M2(b — a)°'/T{a + 1) is independent of |e| < 1. 
Finally, from Inequalities ([78]) . ([80]) and ([8T]) . we have: 

V|e| < l,Vt G [a,b], ||g"+^^'"(t) - g'''"(t)|| < Mr f||g"+^"'" - g"'" 



(t)<M3, (81) 



(t) + M3|e|. (82) 
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Using the fractional Gronwall's Lemma \W[ we conclude that: 

V|5|<1, VtG[a,6], ||g"+^"'"(t)-(?"'"(t)|| < Ms\e\E^,i{Mit-ar) (83) 

< M3|e|E,,i(M(6-a)°). (84) 

Defining Ci := MsEa,i{M(b — a)"), the proof is completed. D 

A. 3. Result of stability of order 2. — In this section, we use the notations and definitions 
given in Section 12.21 We prove the following Lemma with the help of Lemmas [16] and [T71 

Lemma 18. — Let u, u £ 'rf^{[a,b],W"'). There exists a constant C > such that: 

y\e\ < 1, Vt G [a,b], |k"+''''°(t) -g"'"(t) -eg(t)|| < Ce\ (85) 

where q is the unique global solution of the following linearised Cauchy's problem: 



df df 

,D^q=-^{q''''',u,t) xq+-^{q'''»,u,t) xu 

q(a) = 0. 



(LCP^ 



The existence and the uniqueness of q are given by Theorem [7] and Condition ( fx lip ) . 



Proof. — We proceed in the same manner than for Lemma [T71 We have for any |e| < 1 and 

any t G [a, h\: 



|g"+"^'°(i)-g"'°(i)-eg(t)|| 
<r ( /(g-+-'",n + en,t)-/(g"'-,n,t)-e^(g"'°,n,t)xg-e^(g-'°,n,t)x^ 



dx 



dv 



it). 
(86) 



With a Taylor's expansion's of order 2 with explicit remainder, we have for any |e| < 1 and 

any y e [a,b]: 

\\f{q^+^^'-{y),u{y)+eu{y),y)-f{q'^'"{y),n{y),y) 

-e^{q-'"{y),uiy),y)xqiy)-e^{q^'^{y),u{y),y) x u{y) 



< 



^^ {q^'''{y),u{y),y) x (g«+^^'°(?/) - g"'"(y) - eqiy)) 



dx 



+ 



^y^ f{ei{y)My).y){'f^''''%y) - 'f'%y),eu{y),Qy 



(87) 



where for any y G [a, 6]: 

- Cf(y) e [g"'"(y),g"+^^'°(y)] C [|k'''"||oo - Ci, ||g"'"||oo + CiY according to Lemma [HI 
Then, there exists Mi > such that ^f(y) € [—Mi, Mi]"' for any |e| < 1 and any 

y G [a,b]; 

- CKy) G [u{y),u{y) + eu{y)] C [-M2,M2]"' with M2 > independent of |e| < 1 and 
y € [a, b] since u and u are continuous on [a, b]. 
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Since V^/ is continuous on the compact [— Mi, Mi]'^ x [— M2, M2]"* x [a, 6], there exists M3 > 
such that for any |e| < 1 and any y G [a, 6], we have: 



^V2/(ef(y),e|(y),2/)(g"+^^'"(y)-g"'"(y),eu(y),0)^ 



<M3(||g"+=^'°(y)-g"'' 



+ 2\\q 



,M+£M,a 



(y)-g"'"(y)lllk^ 



+ eii 



')■ 



Then, from Lemma [T71 and since u is continuous, there exists Af4 > such that: 

1. 



V|e| <1, VyG [a, 6], 



TjVV(^!(y), ei(y), y) (9"+^^'"(y) - g"'"(y), en(y), o)^ 



< M4e^ 



Consequently, from Inequahty (f87|l and since / satisfies Condition ( /x HpD (implying 
\\df /dx\\oo < Af), we have for any |e| < 1 and any y G [a, 6]: 

||/(g'^+^"'°(y), u(2/) + eu(y), 2/) - /(^"■"(y), u(2/), y) 



. / ( „u,a 



dx 



q'''"{y),u{y),y) x q{y) 






< 2M||g"+-'"(y) - g-'"(y) - eg(y)|| + M^e'. (90) 



Finally, using the previous inequality and Inequality ()86p . we obtain for any |e| < 1 and any 

t e [a,b]: 



u+eu,a 



(t) 



'{t)-eq{t)\\ <2MI 



a ( \\ u+eu,a „'^,oi 



q'''"-eq\\){t)+M,e' 



(91) 



where M5 := M4(6 — a)°/r(l + a) is independent of \e\ < 1. Finally, from the fractional 
Gronwall's Lemma [T6l we conclude that: 

V|e| < 1, Vt G [a,b], |k"+''''"(t) - g"'"(t) - eg(t)|| < Af5e2E„,i(2M(5 - a)"). (92) 

Defining C := M5Eo,i(2M(6 — a)"), the proof is completed. D 

A. 4. Proof of Lemma [5l — In this section, we prove Lemma [5] and consequently, we use 
the notations and definitions given in Section [2.21 



Let u,u£ ^°([a,6],M™) and q G ^^([a, 6],IR'^) the unique global solution of ( |LCP|p . From 
Lemma [THl we have: 



Vkl < 1, q 



,u+eu,a „«,o 



(7"'" + eg + /l^ 



(93) 



where ||/i^||oo ^ Ce^. In particular, since q is continuous, there exists Mi > such that: 

V|e| < 1, Vy G [a, b], [(/"'"(y), g"+^"'"(y)] C [-Mi,Mif. (94) 

In the same way, since u and u are continuous, there exists M2 > such that: 

y\e\ < 1, Vy G [a,b], [u{y) , u{y) + eu{y)] C [-M2,M2]™. (95) 

We have: 



V|e| < 1, C^iu + eu) - £°(n) = / L(g"+^"'°, -u + eu, t) - L(g"'", u, t) dt 



(96) 
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With a Taylor's expansion of order 2 with exphcit remainder, we have for any |e| < 1 and 

any y e [a,h]: 



Liq 



< 



u+eu,a 



{y),u{y) + eu{y),y) - L{q''''^{y),u{y),y) 

dv 



fir 

^1^ ('?"'"(?/), n(y),y)-g-(y) 



f''^iy),u{y),y).u{y) 



dL 



— {q-'-{y),u{y),y)-h^y) + -V^L{ei{y),e2{y),y){eq{y) + h'{y),eu{y),Oy , (97) 

where for any \e\ < 1 and any y G [a,b], ^f(y) G [g"'"(y),g"+^'"'°(y)] C [—Mi, Mi]"' and 
^2(2/) S [■'^(y)i ^(y) + £^(y)] C [— M2, M2]'^. Since L is of class ^^, we obtain easily that there 
exists M3 > such that for any \e\ < 1 and any y S [a, b]: 



Liq 



u+eu,a 



(y), u{y) + £«(?/),?/) - L{q''''^[y), u{y),y) 



dL 



dL 



e^{q''^''{y),u{y),y)-qiy)-e—{q^'^{y),uiy),y)-u{y) < M^e\ (98) 



Consequently, we have for any < |e| < 1: 



Hence: 



lim 



£"(-« + eu)-£"(u) 



(g"'", n, t) • ^ + ^(g'^'", n, t) • -u dt. 



9x 



9f 



< Mge. (99) 
(100) 



The proof is completed. 
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